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Diffuse Gamma Radiation 


C. E. Fichtel, G. A. Simpson*, and D. J. Thompson 
NASA/Goddard Space Flight Center, Greenbelt, MD 20771 

ABSTRACT 

An examination of the intensity, energy spectrum, and spatial 
distribution of the diffuse /-radiation observed by SAS-2 away from 
the galactic plane in the energy range above 35 MeV has shown that it 
consists of two- components . One component is generally correlated with 
galactic latitudes, the atomic hydrogen column density as deduced from 
21 cm measurements, and the continuum radio emission, believed to be 
synchrotron emission. It has an energy spectrum similar to that in the 
plane and joins smoothly to the intense radiation from the plane. It 
is therefore presumed to be of galactic origin. The other component is 
apparently isotropic, at least on a coarse scale, and has a steep energy 
spectrum. No evidence is found for a cosmic ray halo surrounding the 
galaxy in the shape of a sphere or oblate spheroid with galactip 
dimensions. Constraints for a halo model with significantly larger 
dimensions are set on the basis of an upper limit to the y-vay 
anisotropy. - 

The features of the galactic component of the /-radiation and con- 
clusions which can be Inferred from the data are:, (1) A linear correla- 
tion between both the 35 MeV < E < 100 MeV and the > 100 MeV /-ray 
intensity and the radio continuum radiation as measured at 150 MH^* 

(2) For jb| > 12.8°, a linear correlation between the /"ray intensity 
in the same two energy ranges and the deduced atomic hydrogen column 


*NAS-NRC Postdoctoral Research Associate, 
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density. (3) For |b| < 12,8® and 60® < & < 300®, a correlation between 
the y-ray intensity and the atomic hydrogen column density deduced from 
the 21 cm emission, but with less y-radiation than would be expected 
from the linear extrapolation .based on data for- jb-j- ■> 12'.8®. (The 
magnitude of the decrease suggests that both the cosmic ray density 
and the molecular- to- atomic hydrogen ratio decrease beyond the local 
region, < 0,5 kpc, for galactic radii larger than that for the solar 
system, although the possibility that only one of these two components 
decreases dramatically cannot be excluded.) (4) For 300® < -S' < 60®, a 
7 -ray intensity from the galactic plane above that suggested' by a 
linear extrapolation based on data for jbj > 12.8°. (5) Ah energy 

spectrum of the high latitude y-rays which is 'similar to that from the 
central galactic plane and hence has a substantially larger intensity 
in the 35 Me'V to 100 MeV range than that expected from -most current 
estimates of the interstellar cosmic ray electron spectriim. If this 
35 MeV to 100 MeV excess emission is interpreted as being due largely 
to cosmic ray electron bremsstrahlung, a significant enhancement in the 
low energy (35 MeV to several hundred MeV) cosmic ray electron spectrum 
over previously assumed spectra is implied. If the excess is attributed 
to Compton radiation, the electron intensity at > 50 GeV must be sub- 
stantially higher than current measurements at the earth suggest. (No 
enhancement is necessary for the 1 to 30 GeV range, so there is not 
necessarily a conflict with the radio continuum observations . ) 

The "isotropic” component was found to have the following properties; 
(1) an energy spectrum which is quite steep, with a 2.85lo!35 differential 



3 


power law index for a best fit straight line between 35 MeV and ~ 150 

MeV., (2) an intensity above 35 MeV of (6. 3+1. 4) xlO photons 
-2 -1 -1 

cm' ster s , and (3) an extrapolated intensity at 10 MeV which agrees 
well with the diffuse ."isotropic" intensity measured at low energies. 
Relatively few theoretical models are consistent with the combination 
of the observed spectral shape, the intensity, and the lack of strong 
anisotropy, and those are not entirely free of other difficulties. 

Key words: galaxies: Milky Way - galaxies: structure - 

gamma rays : general 

I . INTRODUCTION 

There is substantial interest in the diffuse y-radiation from 
regions away from the galactic plane because of the several possible 
origins of this radiation, some involving important cosmological 
implications and others relating to the nature of our Galaxy. The 
possible explanations of this radiation include emission from normal 
or extraordinary galaxies of various types, interactions of the universal 
black body radiation and cosmic rays (either primordial or from galactic 
leakage), primordial black hole emission, particle-antiparticle inter- 
actions at the boundaries of superclusters of matter and antimatter 
resulting from the baryon symmetric big bang theory, as well as various 
galactic models involving primarily brerasstrahlung and Compton radiation. 
These latter theories would in general suggest a variation of intensity 
with galactic latitude, 

A large number of experiments have been conducted in an attempt 
to study the diffuse radiation. At low y-ray energies, the situation 
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is rather unclear, in part due to the background problem which has been 
encountered in this energy' range. At intermediate energies (from about 
8 to 35 MeV) , the most serious difficulty is probably the lack thus far 
.of an .adequate, .exposure .opportunity outside the -Earth* s .atmosphere . 

At higher energies (E > 35 MeV) , which is the range to be considered 
here, the Second Small Astronomy Satellite (SAS-2-)' has now provided, 
data which confim the existence of a diffuse high latitude y-radiation 
(Fichtel, Kniffen, and Hartman 1973) first suggested by the OSO-3 
7 ~ray data (Kfaushaar et al. 1972). 

In a letter summarizing the initial analysis of the entire SAS-2 
data set Fichtel et al. (1977) showed that when the existing SAS-2 high 
energy data were stumnarized over galactic longitude, there was a strong 
variation with galactic latitude. For the region 10° ^ 1^1 ^ 90°, it 
was shown that the summed data were consistent with an expression of 
the form Cj_ + C 2 /sin |b|, corresponding to an isotropic component plus 
a galactic disk component. The ratio was larger at 35-100 MeV 

than at energies above 100 MeV indicating that the apparently isotropic 
component on a coarse scale had a steeper energy spectrum than the 
galactic component. The galactic component Intensity joined smoothly 
to the more intense region along the plane reported earlier (Fichtel 
et al, , 1975) and had a similar energy spectrum. 

In this paper, -a substantially more detailed analysis of the 
diffuse radiation is performed along with a comparison to relevant 
data at other wavelengths, including the 21 cm emission from atomic 
hydrogen, the radio continuum radiation (a large portion of which is 
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believed to be synchrotron radiation) , and the limited ultraviolet and 
radio information on .the local molecular hydrogen. It is found that 
there are indeed two quite distinct components to the diffuse radia- 
tion, one which shows a glood correlation with the galactic matter 
distribution and the synchrotron radiation^ and the other, having a 
much steeper energy spectrum, which appears to be isotropic at least 
on the scale that it could be examined here. The galactic component 
is interpreted in terms of its implications for both the local and 
more distant regions of the galaxy. The apparently isotropic 7 -radiation 
is then discussed particularly with regard to the constraints placed on 
possible models by the steep energy spectrum, the intensity, and the upper 
limit on the anisotropy. 

II, OBSERVATION AND ANALYSIS 

The principal data sources are: (1) the 7-ray s"ky survey of SAS-2, 

in the energy intervals 35 < E < 100 MeV and E > 100 MeV to he dis- 
cussed in this paper, (2) the 21 cm summaries of Daltabuilt and Meyer 
(1972), Heiles (1975), and Burton (1976) which provide a measure of the 
atomic hydrogen column density, and (3) the 150 MH^ synchrotron radia- 
tion survey of Landecker and Wieleb inski (1970). In addition, to esti- 
mate the molecular hydrogen component, the results of the Copernicus 
ultraviolet survey of Savage et al. (1977) and the recent radio survey 
of the 'CO emission line analysis by Gordon and Burton (1976) are used. 

The SAS-2 7-ray telescope is a 32- level thin plate digitized spark 
chamber o A detailed' description of the instrument is given by Derdeyn 
et al. (1972). The entire SAS-2 celestial 7-ray data set, which was 
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collected during the period from November 1972 through June 1973, 
has now been analyzed in accordance with procedures described by 
Fichtel et alo (1975) « The analysis used the detailed sensitivity, 
angular response function, y-ray arrival direction accuracy deter- 
mination, energy, measufeinents., and energy resolution, function (for 
determining the true energy spectral shape) determined in the 

extensive calibration outlined in that paper, as well as extensive 
(> 80%) rescans of the y-ray event films to search for possible in- 
efficiencies and selected earth albedo measurements during the satellite's 
history to check for possible changes in detector performance. Figure 1 
shows the region of the sky covered by the complete SAS-2 y-ray observa- 
tions. Notice that a full strip between the plane and the north pole 
exists in the vicinity of = 0° and also for H « 240°, or in the general 
anticenter region well away from the intense central region of the 
galactic plane. 

In the analysis to.be presented here, the galaxy has been divided 
into eight parts, which are shown in Figure 1. Region A encompasses 
the inner galaxy, regions B and F allow study of directions tangent to 
the solar circle, and regions C through E display the anticenter. Each 
of regions A through E includes all latitudes with jb] < 60° and all 
longitudes for which SAS.-2 y-ray observations exist with the exception 
of the regions in which the strong contributions from the Vela pulsar 
PSR 0833-45 (Thompson et al. 1975 and Thompson et al. 1977a), Cygnus 
X-3 (Lamb et al, 1977b), .the Crab pulsar PSR 0531+21 (Kniffen et al, 

1974 and Thompson et al. 1977) and the y-ray source at ^=195°, b=+5° 
(Thompson et al. 1977b) exist. Regions for which b > 60“ and b < -60“ 
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(G and H respectively) have also been included in our analysis. The 
choice of the analysis regions -was determined to a significant degree 
by the availability of 7 -ray observations and the photon statistics. 

For purposes of orientation, the E > 100 MeV and the 35 MeV < E 

< 100 MeV 7 ~ray intensities are plotted as a function of galactic 
latitude for regions A through F in Figures 2 and 3. (See the caption 
to Figure 2 regarding the reason for the latitude interval boundaries 
not being integers.) Notice first that there is a significant variation 
with latitude. 

Note also in Figure 2 that near J0=O°, there is an excess on the 
positive latitude side of the plane in the 60° < £, < 100 ° and the 
100° < I < 150° intervals and on the negative side in the 210° < I 

< 250° interval. This effect is similar to the "hat brim" effect seen 
in the galactic matter distribution, and is presumably a reflection 
of large-scale galactic structure. 

In order that they may be compared to the 7 -ray data, the column 
densities of atomic hydrogen, and the 150 MH^ brightness tempera- 

ture, l]_ 3 Q,have been averaged over longitude within narrow latitude 
intervals for each of the regions in which 7 -ray data exist. For . 
direct comparison to the 7 -ray data, these intensity measurements have 
been "defocused" to take into account the lesser angular resolution of 
the 7 -ray observations. The knowledge of the molecular hydrogen in the 
local region is not sufficient to make a similar comparison, and it will 
be discussed later i 

The correlations between the atomic hydrogen column density and 
both the low energy and high energy 7 -ray emission are shown in Figures 
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4 and 5 for jb| > 12,8°, In. these figures, points near the origin, 
having small column densities, correspond to high latitude j-ray obser- 
vations, since the shortest paths out of the galactic disk are at high 
latitudes. Eor the .same, reason, these points represent predominantly 
local matter. The low- latitude data correspond to high column densities, 
and include radiation from the more distant regions in the survey. 

The local y-ray matter correlation is described rather well by a 
linear relationship of the form: 

= ANiji + ^ Cl) 


where the column density of atomic hydrogen in units of atoms/cm^, 

and is the /-ray intensity in units of protons cm""^ster" ’^s”'*'. ^0 

straight line shown on each sub- figure of Figures 4 and 5 is the least 
squares fit to the combined observations from all the sby intervals 
with |b I s 12.8°. Each fit has \ = 0.8 where v is the number of 

degrees of freedom. In view of the thickness of the matter disk, the 
latitude of 12,8“ corresponds to a distance of about 500 pc- to the 
point of one scale height from the plane. Within the uncertainties 
of the measurements all of the regions are seen to be consistent with 


this single curve. The values of "A" and ”B" in equation (1) are: 

“26 “1 
3o0 X 10 photons cm“^ster"ls"l(atoms/cm^O and 1.0 x 10“-^ photons 

cm”2stej^“ls-l for E > 100 MeV and 4.3 x 10 photons cm”^ster“ls“l 
(atoms/cm^) gjid 4.9 x 10 photons cm~^ster ^s ^ for 35 ^eV < E < 


100 MeV. The uncertainties in these numbers are illustrated in Figures 
4 and, 5 since the coupling of the errors precludes giving a single 
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xneaningfui one-parameter error estimate. Clearly the first term 
represents a source which is proportional to the matter density and 
the second term represents an isotropic component. The 35 MeV < E 
< 100 MeV interval has a relatively larger ”B" term, implying a steep 
energy spectrum for the isotropic component, at least on a coarse 
scale . 

The low latitude intensities, involving integrals over more distant 
points, show significant deviations from the linear relationship. In all 
regions for which data exist free of substantial clearly established point 

source contributions except A (See Pigure 1) , the low latitude points 
fall below the previously determined line, as seen in Figures 6 and 7. 

In these regions, the y-ray intensity from the low latitudes corresponds 
to a substantial contribution from regions lying at distances from the 
galactic center greater than that of the Sun. By contrast, in region A 
the y-ray intensity within a few degrees of b=0° lies above an extrapo- 
lation of the same straight line. 

For comparison to the continuum radiation, the 150 MH^ survey 
of Landecker and Wielebinski (1970) was used. The 150 MH 2 radiation 
has characteristics which are very similar to those of the y-ray 
emission. Figure 8 shows the 150 MH 2 observations plotted as a function 
of the 21 cm column density. For jb] s 12.8° the relationship between 
the brightness temperature and the matter density may be well repre- 
sented(x^/v = 1-2) by = (1. 55x10” +165. The uncertainties 

in these numbers are given in the figure. Since the temperature 
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is directly proportional to intensity, this expression is identical in 
form to that of equation -(1). for regions away from the galactic center, 
the low- latitude 150 MHL^ data points fall below the line, in the same 
manner as the y-radiation Illustrated in Figure 6= 

The correlation between the 150 continuum radiation and the 
35-100 MeV y-ray data is shown in Figure 9. This figure shows a good 
linear correlation between the y-ray and 150 intensities 
with no evidence of a significant deviation from linearity or a non- 
zero intercept. The least squares fit' straight line and its error are 
shown in the figure. The E > 100 MeV y-ray intensities show a similar 
correlation. In addition, on a large scale the y-radiation, the 21 cm 
column density, and the 150 MH 2 brightness temperature are all reasonably 
consistent with a function of the form C]_+C 2 /sin jb} for ]b ] > 10°, as 
noted previously for the y-rays (Fichtel et al. 1977). and for the 21 cm 
column density (e.g. Heiles, 1976). This relationship is in fact to be 
expected for any form of diffuse radiation originating in the disk of 
the galaxy, simply from the geometry of a thin disk. 

Finally, the energy .distribution of the diffuse radiation was 
examined as a function of latitude and longitude. These data are shown 
in Figure 10. For this comparison it was necessary to multiply the high 
energy (> 100 MeV) y-ray intensities by an appropriate resolution 
function (thereby in effect defocusing the high energy data) so that 
they would be directly comparable to the 35-100 MeV y-ray measurements, 
for which there is a poorer angular resolution. A least squares analysis 
of these data shows that a linear correlation is consistent with the 
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data, as shown in Pigure 10, (x^/v=0o8). The straight line, in fact, 
fits the data within uncertainties over the entire latitude and longitude 
range included in the data set. There is no evidence for a variation in 
7 ~ray energy spectrum with direction, although within present uncer- 
tainties the possibility of more than one major source with a somewhat 
different energy spectrum cannot be eliminated, Finally, the intercept 
not being at (0,0) is another indication of the different spectral shape 
of the isotropic component. 

Because of the limited energy resolution of the detector the 
energy spectrum cannot be determined in fine detail. Further, the 
previously determined constants, ”A” and for equation (1) cannot 
be used directly in determining the constants of an assumed spectral 
shape, but must be used together with the energy resolution function 
(Fichtel et alo, 1975) and an assumed spectrum. If power law spectra 
are assumed, for the "isotropic" component a differential index of 
2.85^8ri§ and intensities above 35 MeV and 100 FfeV of (6.3+1.4)xl0“^ 
and (0.9+0.4)xl0"^ photons cm~^s~^ster"^ respectively are obtained. 

For the galactic component, a differential spectral index of lo7H^o3 
and intensities above 35 and 100 MeV of (6 , 9+1.7)xl0”^^ and 

(3 .3i4) . 8)xl0 ^ster ^ respectively are obtained. 

-5 5 

These correspond to intensities of (1.6 +0.4)xl0 and (0.8;Wo2)xl0 
- 2 - 1-1 

photons cm s ster for galactic latitudes near the pole and 
(8.3+2,l)xl0 and (4o0 +1.0)xl0 ^ photons cm''^s~^ster”^ for a typical 
region with a galactic latitude of about 15® » l£ 3 combination of a 
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a nuclear interaction and power law spectriJin (of the Compton or 

bremsstrahlung type spectral index) is assumed, almost identical results 
are obtained primarily because a relatively large power law component 
must be assumed. The two deduced spectra are shb\ra individually together 
with their sum' in Tigure 11. These two spectra are very similar (agreeing 
within the indicated uncertainties) to those determined earlier assuming 
^1 C2/sinjb| dependence for the y-rays (Fichtel et al. 1977). 

Finally, there is the matter of calculating the degree of isotropy 
of the ’’isotropic" component. To accomplish this goal, the galactic disk 
component has to be subtracted first in a manner which is as unbiased as 
possible. The material already presented has shown that the atomic 
hydrogen colxmin density, the synchrotron radiation, and the y-radiation 
are all linearly correlated with each other for Jb| > 12,8° on the coarse 
scale being considered here. In the next section, it will be seen that 
because of the nature of the distribution of the two components involved 
in each of the interactions believed to be relevant, theoretical con- 
siderations also indicate that the major components of the local galactic 
y-radiation should be proportional to either the matter column density 
(high energy, nuclear interactions, and bremsstrahlung) or the cosmic 
ray electron column density (Compton radiation) which is approximately 
represented by the continuum radiation. Since the matter column 
density has no known component beyond the galactic disk whereas the 
continuum radiation may very well have, the previously determined linear 
correlation between the HI column density and the y-ray intensity was 
used to determine the multiplying factor for the colttmn density. 
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The galactic disk 7 -ray component determined in this manner was then 
subtracted from the total to find the 7 -ray component not associated 
with the galactic disk in the various regions of the sky being con- 
sidered „ This approach does attribute a small center-to-anticenter 
dependence to the' galactic disk 7 -ray component, incidentally. 

The ratio of the average residual (and larger portion) of 
the 35-100 MeV 7 -ray intensity in the directions (300° < A < 60°, 
20° < Ibl < 40°) to that from (100° < H < 250°, 20° < |bj < 40°) was 

then determined. The ratio was found to be 1.10 dp. 19* with a cor- 
*For the reader who may be interested, the similar ratio for the total 

Y~ray intensity before subtracting the galactic component is 1.24+0.15 

responding 987;, confidence upper limit of 1.48. Hence, there is no 
evidence for a residual center-to-anticenter anisotropy. The ratio 
between the average 7 -ray intensity from regions with jb ] > 60° to 
that from 20°<|bj < 40° is found to be 0.87 +0.09. This result 
naturally argues against a thin oblate spheroid. Both of these ratios 
will be discussed in terms of constraints on specific models in 
section III. 

IIIo DISCUSSION 

From the material presented in the previous section it is 
apparent, both on the basis of the latitude and the energy distribu- 
tions, that the 7 -radiation consists of two components, one galactic 
based on its latitude and longitude distribution and its correlation 
with the matter distribution, and the other apparently isotropic, at 
least on a coarse scale, with a relatively steep spectrum. These two 
components will nox^ be treated separately. 



a) The Local Galactic Component 
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Based on its angular distribution and intensity, the y-ray . emission 
from the galactic plane has generally been assumed to result primarily 
from cosmic ray interactions with matter, with a much smaller contribution 
from localized sources. The Tatter are discussed' lafer'Th fh'is subsection. 
Calculations related to other production mechanisms in the local galactic 
region including Compton radiation from cosmic ray electrons interacting 
with starlight, infrared photons, and the blackbody radiation show that 
these may also make a contribution, but that synchrotron radiation is 
quite negligible (e.g, Schlickeiser and Thielheim, 1976; Worrall and Strong, 

1977; Fichtel et al. , 1976; Piccinotti and Bignami, 1976; and Stecker, 1977). 

It seems reasonable to test the assumption that the latitude dependent 
portion of the high-latitude y-ray emission also has its origin principally 
in cosmic ray matter interactions in view of the similar energy spectrum, its cor- 
relation with the 21 cm emission, and its smooth connection to the 

radiation from the galactic plane. This hypothesis, which will now be 

pursued, will be seen to have some difficulties with the observed 

intensity and possibly the energy spectrum. Explanations for the 
differences will be suggested. 

^ line of sight in the direction (j2t,b) for the range of para— 
ms ter s being considered here, the y-ray intensity produced by interactions 
between cosmic rays and matter can be shown (Fichtel et al, 1976) to be 
given by the expression: 

I(E^,^,b) = dr [Sy^(E^,r=0)g^Cr,je,o)£^(r,£,b) 

+ (E ,r=0)gg(r,£,b) fjj^(r,X,b) (2) 

© ' 

P 

+ (E^,r=0)g^(r,^,b) f^^(r,.6,b)] 
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where f is the ratio of the total interstellar gas density at a- dis- 


tance r from the sun in the direction (j&jb) to -that at r=0'. S' repre- 

^n 

represents the y-rays produced per second per cm in interactions of 
nucleonic cosmic rays (with the intensity and spectral distribution in. 
the solar vicinity) with the interstellar gas, and S and S are 


p s 

similar functions for primary and secondary cosmic ray electrons, 
respectively, and express the spatial variation with galactic 
position of the ratio of the primary cosmic-ray nucleon and electron 
components respectively to their interstellar value in the solar 
vicinity. 

It is seen readily from equation (2) that, if the cosmic ray 
primary electron component is proportional to the cosmic ray nucleon 
component, as is commonly assumed, and if the secondary electron 
component is small compared to the primary electron component locally 
as the positron data suggest, the y-ray intensity is simply proportional 
to the integral over the product of the ratios of the cosmic ray inten- 
sity and matter density to their local values. 

Hence equation (2) may then be rewritten in the approximate form 
I(E ,X,b) p dr S (E^,r=0) g (r,^,b) fjjj(r,^,b) (3) 

where the subscript cr refers to the combination of the cosmic ray 
electrons and protons. 

With regard to the matter, as has already been noted, a reasonably 
accurate estimate of the atomic hydrogen column density in any direction 
can be obtained from the 21 cm measurements; however, the current status 
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of the knowledge of the other major component, namely the molecular 
hydrogen, is much less- clear. A detailed picture of the molecular 
hydrogen density or even the column density does not exist. From the 
Cbpernicus data on the local interstellar gas densities, Savage et alo 
(1977) estimate the fraction of the gas in molecular form to be at 
least 0.25 and possibly as high as 0.5. Although not directly applicable 
to the local region, the large-scale survey of Gordon and Burton (1976) 
indicates that the fraction of gas in molecular form may be about 0.5 
for an average distance of 10 kpc from the galactic center. Thus a. 
reasonable range for the average ratio of the molecular hydrogen density 
to the atomic hydrogen density in the local region of the galaxy appears 
to be 0.3 to 1.0. There is the additional consideration that molecular 
hydrogen .is generally believed to be preferentially concentrated in 
clouds; however, if, as will be discussed later it is true that the 
cosmic ray density is uniform on a sacle of the average cloud- to-cloud 
dimension, this difference in distribution between the atomic and 
molecular hydrogen is not of concern here. Other tracers of the total 
matter, such as interstellar reddening and galaxy counts, are 

generally correlated with the HI column density in a linear manner for 
the local region (Heiles, 1976), although there are some small scale 
differences, 

A consideration in relation to the bremsstrahlung radiation is 
the percentage of helitam nuclei and heavier nuclei in the interstellar 
medium, because of the Z (Z-fl) dependence of this radiation. Fichtel 
et al, (1976) assumed the heliura-to-hydrogen ratio to be 0.10 and the 

i 
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heavy- to -hydrogen ratio to be 0.01 leading to contributions relative 
to hydrogen of 0.30 and 0.25 respectively. There seems to be no basis 
for changing these estimates, although there is substantial uncertainty 
in them. If the interstellar medium were particularly rich in heavy 
nuclei, the bremsstrahlung could be somewhat larger, but in fact a 
heavy-to -hydrogen ratio of 0.01 appears to be on the high side of 
current estimates. 


Turning now to the y-ray source functions, in equation (3) 

cr 

consists of the sum of and , where in turn is the sum of , 

' ^ -a. '^e 

P 

the primary electron contribution, and a small correction for the 

secondary electron term. For S (E > 100 MeV) there is the recent 

^n 

” Z6 

value of 13,0x10 photons (E > 100 MeV) cm'^s"^ from the work by 
Fichtel et al. (1977b) and Kniffen, Fichte 1, and Thompson (1977) based 
on the earlier calculations of Stecker (1973), but using more recent 


values for the parameters. For (E > 100 MeV) the same authors give 

e 

— 26 "*3 “X 

3,5x10 photons cm s . These numbers are based on there being 1«04 
protons and electrons, or ’’atoms”, per cm"^ in either atomic or molecular 


form, a helitun-to-hydrogen ratio of 0.1, and a heavy-nuclei-to-hydrogen 

“ 26 

ratio of OoOl. Using the same procedures, values of 3.3x10 photons 

cm ^s ^for S (35 MeV < E < 100 MeV) and 5.9x10”^^ photons cm for 

(35 MeV < E < 100)MeV are obtained. More recently Savage et al. 

® “3 

(1977) have estimated the number of ’’atoms” cm to be 1.15 or possibly 

slightly larger, and hence (E > 100 MeV) and S (E > 100 MeV) would 
-26 ^ 

become 14.4x10 and 3.9x10"^® photons respectively, and 
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S (35 MeV < E < 100 MeV) and S . (35 MeV < E < 100 MeV) become • 

e 

3.6x10 ^^photons cm ^ and 6 „5xl0”^®pliotons respectively. 

Hence, (E > 100 MeV) becomes 18,3x10 ^^'photons cm ^ and 
'cr 


Sy (35 MeV < E < 100 MeV) becomes 10.1x10”^^ photons cm“^s"'*‘. If it 
cr 

is assumed that there might be as much as a factor of three more 

molecular hydrogen, i,e, that the molecular-to-atomic hydrogen ratio 

is 1.0 rather than 0.3, then S (E > 100 MeV) and S„ (35 MeV < E 

cr ' cr . ■ 

< 100 MeV) might be as large as 27.5xl0"2^ photons cm'^s"! and 

15.2 X 10 photons cm~‘^s”''' respectively. 

In addition to the uncertainty in S due to the uncertainty in 

^cr 

the molecular hydrogen density, there is an uncertainty in the electron 

bremsstrahlung contribution (which results primarily from electrons in 

the 10 to a few hundred MeV range^ which in turn is determined largely by 

the lack of knowledge of the electron cosmic ray spectrum in the apprapriate 

energy range. From an analysis of these uncertainties, Shukla and 

Cesarsky (1977) conclude that an upper limit to the electron source 

function is about 50% larger than the source function used here. Such 

a picture would be consistent with the closed galaxy model (Rasmussen 

and Peters, 1976). On the other hand, Webber (1977) has argued on the 

basis of radio data that the interstellar electron spectrum (in the 

energy range above 1 Ge'^ is actually lower than the stardard demodulated spectrum 

(Cummings, 1973). The y~ray spectrum observed here and that observed in the 

galactic plane (Fichtel et al. 1975; and Bennett et al. 1977) are consistent 

with a larger electron contribution or some other similar power law spectral 

component. If the possibility of a 50% larger electron intensity is 


,“ 3-1 
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considered j S (E > 100 MeV) and S (35 MeV < E < 100 IfeV) would have 

“26 "3 “1 “26 

values of 30.4x10 photons cm and 20.1x10 photons cm"^s~^ 

respectively. Hence, from these considerations, the range of values 

for (E > 100 MeV) and (35 MeV < E < 100 MeV) are (1.8 to 3,0) x 
'^cr '^tr- 

io photon cm ^ and (1.0 to 2,0)xl0~^^ photon cm~^s”^respectively. 

If the cosmic ray density in the local region is essentially 
constant near the plane (The justification for this assumption will he 
discussed in the next paragraph.) where the great majority of the cosmic 

ray matter interactions occur, then g of equation (3) becomes 1,0 and 

cr 


equation (3) becomes: 


7 . 




(4) 


Assuming that in the local region the ratio of molecular hydrogen to 
the total matter is invariant on a broad scale averaged over a column 
as discussed before, equation (4) becomes 


4tt P 


(5) 


where Njj^(X,b) = J* (r,J&,b) dr is the column density of atomic 
hydrogen which can be determined from 21 cm radio observations and P 

s 

is the local atomic hydrogen density. p is taken to be ,0.86 atoms 

O 

cm to be consistent with the value used in determining . Then 

^cr 

. -26 
I (E > 100 MeV,J0,b) = (1.7 to 2.8)xl0 N (j6,b) 


7 


cr 


HI 


( 6 ) 
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" 26 

I (35 MeV < E < 100 MeV,J6,b) = (0.9 to 1.9)xlo" (7) 

^cr 

There are several considerations supporting the point of view that ' 
the co,sm'ic ray distribution is relatively uniform in the' vicinity 
(^0.5 kpc) of the solar system. First, the observed cosmic ray 
anisotropy is extremely small. Further, whereas the cosmic rays are 
expected to be correlated with matter on the scale of galactic arras 
within the plane because the cosmic rays and magnetic fields on a 
broad scale can only be held by the gravitational attraction of the 
matter (Parker, 1966 and 1969), the cosmic rays are expected to be 
uniform on a cloud- to- cloud scale within the plane, although varying 
with height above the plane (Kniffen, Fichtel, and Thompson, 1977). 

The concept of local uniform cosmic ray density is also supported by 
the analysis of Freier, Gilman, and Waddington (1977) who show that the 
level of the y-ray emission can only be explained if the cosmic rays' 
are uniform locally since the predicted y-ray emission woul'd be much 
too large for the cosmic rays to be proportional to the- matter on tbe 
scale of clouds. Further, since the magnetic fields are expected to 
bulge significantly between the clouds, the scale height of the cosmic 
rays is expected to be large compared to that of the matter (e.g. 

Parker, 1966; Bignami et al. 1975; Eniffen et al. 1977). This concept 
is supported by the large scale height (~ 0.75 kpc) observed for the 
radio continuum radiation (e.g. Baldwin, 1976). If it is assumed that’ 
the electron density and magnetic field energy density both have the 
same distribution on the average and that it is gaussian, then since 
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the synchrotron radiation is proportional to the product o£ the two, 
the scale height of each individually is /"2 x; 0g75 hpc» 

This scale height becomes important in the consideration of the 
Compton radiation. Following the work of Fichtel et al, (1977), the 
calculations on which the work was based, and assigning a gaussian 
electron distribution perpendicular to the plane with a scale height 
of /^2 X 0o75 kpc, the Compton radiation resulting from cosmic ray 
electrons interacting with the blackbody radiation is given by the 
following equations: 

^ ^00 MeV) = 0,074 x 10 [sin]b| ] ^ photons cm ^ster ^ (8) 

(35 MeV<E< 100 MeV)=0, 117x10 ^[sinjb j ]“^photons cm''^ster ^s ^ (9) 

bb 

The combined contribution from the infrared radiation and starlight is 
much more difficult to estimate not only because the average photon 
density in the plane is not very well known, but also its variation 
with height above the plane is poorly known. The effective scale height 
appears to be much smaller (~ 0.2 kpc based on the paper of Lambrecht, 
1965) than that estimated for the electrons, so that even though the 
combined source function in the plane is slightly larger (Fichtel et 
al. 1976 and Ficcinotti and Bignami, 1976), the net contribution is 
apparently smaller and is estimated to be: 

1(1 _ (E >100 MeV)=0. 025x10"^ (sin] b j ) ^photons cm'^ster~^s“^ (10) 
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(35 MeV< E< 100 MeV)=0o039xl0 ^(sin|bl)"^ photons cm“^ster 
,ir ... 

Hence, on this basis the combined estin^te for Compton radiation .is.: 

I (E >,100 MeV) = Oo 10x10 . (sin|b|)”^ photons cm ster s (12) 

5 1 

I (35 MeV< E< 100 MeV)=0. 16x10 (sin[bj) photons cm”^ster“^s’^ (13) 

The uncertainty in these last two numbers is quite large, in part due 
to the substantial uncertainty in the electron spectrum in the relevant 
range for the bl'ackbody Compton radiation, namely one to several hundred 
GeV (e.g. Hartman, Ifuller, and Prince, 1977), partially due to the lack 
of knowledge of the electron density distribution perpendicular to the 
plane, and in part due to the uncertainty in the starlight and infrared 
photon distributions o A factor of two uncertainty is tentatively 
assigned to these estimates. 

It is perhaps worth mentioning that the starlight and infrared 

contributions above the plane may be relatively more important in the 

central region of the galaxy, and this possibility should be kept in 

mind when examining low latitudes in the central galactic region. The 

Compton component will be mentioned again in the next section where 

the "isotropic” component is discussed. 

The discussion of the predictions of the intensities is summarized 

in Table I and compared to .the observed intensities using typical N„_ 

HX 

values for a given latitude. In both energy bands, the ranges of 
values allowed by the observations overlaps the upper end of the range 



23 


of values suggested by the theoretical considerations given here. The 
distribution between the two energy ranges as well as the intensity 
does, however, suggest that there is a relatively important and perhaps 
larger Compton or bremsstrahlung type , spectrum contribution since the 
cosmic ray nucleon, matter interaction contribution is relatively 
better knovra and has a spectrum which by itself is inconsistent with 
the observed spectrum. It is important to note that a similar situa- 
tion for both the absolute intensity of the (35 MeV < E < 100 MeV) 
radiation and its ratio to the (E > 100 cfeV) radiation had been, found 
previously for the galactic plane radiation (Fichtel et al., 1975 and 
Bennett et al., 1977). Therefore, this situation is not peculiar to 
the high latitude galactic radiation, but is a general feature of the 
galactic emission at all latitudes. 

In considering the possible origin of this larger than predicted 
galactic radiation (primarily in the lower energy region) three possi- 
bilities seem most likely. These are a large bremsstrahlung contribu- 
tion, a larger Compton component, and a significant point source 
contribution. The latter is a possibility which- is difficult to verify 
or deny with the limited angular resolution of current experiments. 

It does, however, seem a bit ad hoc to postulate a class of low- luminosity 
7 -ray sources with relatively steep spectra to explain these observa- 
tions. The SAS-2 results show no evidence of high-latitude y-ray sources, 
with a typical 2a upper limit of about 2x10"^ photons (E >100 MeV)cm“^s“^, 
and COS-B has also reported no source with a latitude greater than 
|b| =7° (Hermsen et al. 1977). 
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Regarding the possibility of a larger bremsstrahlung contribution, 
it would be possible to have a larger electron intensity in the lower 
energy (35 MeV to several hundred MeV) region relevant to the 7-rays 
being discussed here without demanding a larger intensity in the higher 
energy region which is relevant for synchrotron radiation. It should 
be realized, however, that the required bremsstrahlung spectrum implies 
an electron intensity which is about three to five times that assumed 
by Fichtel et al. (1976) based on the woric of Daugherty et al, (1975) 
in the lower energy region, if it is to be the explanation in itself. 

It is about three if the molecular-to-atomic hydrogen mass ratio of 
1,0 is assumed and the larger number if the molecular-to-atomic 
hydrogen mass ratio is assumed to be one- third. 

If the Compton component is to be the primary contributor, from 
the earlier discussion it is reasonable to assume that. the excess is due 
to a larger than expected cosmic ray electron, blackbody interaction and 
not electron, starlight and infrared interactions since the latter is the 
smaller component and is already defined within relatively narrow limits. 

It is true that the electron spectrum is relatively poorly known above 
about 50 GeV, but there are some reasons for believing that the current 
best estimate used here cannot be exceeded by much. First, the most recent 
spectral measurement (referenced earlier) of many which have been made lies 
a bit below this estimate. Secondly, theoretical considerations on galactic 
trapping and escape are very uncertain, but, if there is an energy-dependent 
escape mode it more likely would lead to relatively fewer high energy 
electrons at greater distances above the plane. There is little 
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experimental evidence to argue against a very thick electron disk of 
lower electron density (thereby avoiding a major conflict with the 
radio synchrotron data) which would increase the volume for the electrons 
and hence the Compton radiation. 

The apparent correlation between the cosmic rays and the matter 
distribution in the spiral arms mentioned earlier might be thought to 
argue against a scale height much larger than the arm dimensions (~ 1 
kpc ) 0 It should, however, be kept in mind that this correlation is 
observed in the inner galaxy, and, in the less defined outer part of 
the galaxy where the cosmic rays may also not be tied to the plane as 
strongly, the scale height for electrons might be larger. 

In view of these considerations, it is worth examining the non- 
thermal radio emission in more detail. As noted earlier, it is generally 
associated with the galactic disk, but has a scale height of about 0.75 
kpc rather than the much smaller ones for galactic atomic and molecular 
hydrogen, .12 and .05 kpc, respectively. The scale height for the . 
continuum radiation is consistent with the concept of the cosmic rays 
and magnetic fields being tied to the plane by the material in the 
more dense regions of the clouds. Clearly the correlation between the 
Compton radiation and the continuum emission would be expected to be 
quite good even though the scale height of the former would be expected 
to be slightly larger as discussed before. 

The physical connection between the continuum radiation and column 
density of matter at high latitudes is not direct, although, as shown 
in Figure 8, observationally there is a strong correlation between the 
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two. The effect may be largely one of geometry, in which roughly disk- 
like distributions of matter and high energy electrons show a correla- 
tion- -because- t-hey both have- -approx-imate-ly -(-sin- b-) latitude dependences 
even though they occupy different regimes in height above the plane. 

This argtiment may also explain the continuum radiation's correlation 
with the disk component of the 7 -radiation (Figure 9), if it is accepted 
that locally the cosmic rays are uniform over the region the matter 
occupies. 

The fact that the best least squares line passes nearly through 
the (0,0) point in Figure 9 is probably fortuitous because both the 
35 MeV < E < 100 MeV 7 -ray emission and the continuum radio emission*' 
have major "isotropic" components which are likely to be of different 
origin, as will be discussed in section IIIc. 

Figure 9 gives a comparison between the intermediate energy 7 -ray 
emission and the continuum emission, which shows that on a coarse scale 
the correlation is in fact generally quite good. The approximately 
linear relationship is to be expected for the reasons outlined in the 
previous paragraph, and the consideration that, although the 7 -ray 
bremsstrahlung radiation comes primarily from electrons in a range 
which barely overlaps in energy, generally falling below that of the 
synchrotron radiation, which in turn falls below that of the Compton 
blackbody radiation, the electron spectrum probably does not vary much 
at least in the local region of the galaxy. 
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b) More Distant Galactic Plane Regions 
Previous studies of the high energy y-rays coming from the galactic 
plane have suggested that either both the cosmic ray intensity and the 
molecular- to -atomic hydrogen ratio increase toward the central part of 
the galaxy or one of the two increases dramatically (e.go Bignami and 
Fichtel, 1974; Paul et al, 1974; Schlickeiser and Thielheim, 1974; 
Bignami et al, 1975; Paul et al. 1975; Stacker et al. 1975; Stacker, 
1976; Puget et al, 1976; Paul et al. 1976), The concept that both 
Increase is the more likely assuming the galactic dynamic balance 
described by Parker (1966 and 1969) and applied to the galactic y-rays 
by Bignami and Fichtel (1974) is essentially correct. It has also 
been noted that the intensity of the high energy /-radiation from the 
plane in the region away from the center (270° < ^ < 90°) suggests the 
cosmic ray intensity is less than the local value (Dodds, Strong, and 
Wolfendale, 1975) o With the data presented in Section II, it is now 
possible to look at these considerations and especially the latter in 
more detail. It should be kept in mind, however that unresolved point 

sources may be making a contribution. 

For small latitudes, marked' deviations from a linear relationship 

between the y-ray emission and the 21 cm column density were observed 
in Figures 6 and 7 for those regions where data exist. In all the 
longitude intervals except 300° < j£ < 60°, the y-ray intensity falls 
below the line derived from the high latitude regions. In each of 
these intervals, the low latitude data include a substantial contribu- 
tion from regions which lie beyond the solar circle. These more distant 
regions are generally in the lower density outer galaxy. The observed 
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depression, of the 7 -ray intensity below the line observed for these 
regions is then consistent with the concept that the cosmic ray density 
is coupled to the total matter density on a large scale, and is hence 
lower in these regions. 

On the other hand, in the region encompassing the galactic center, 
the 7 -ray intensity actually lies above the line deduced for the local 
region for low galactic latitudes. Again, this is consistent with the 
concept of cosmic-ray and matter coupling. 

In addition, studies of the 2.6 mm CO line (e.g. Scoville and 
Solomon, 1975; Burton et al. 1975) indicate that the molecular hydrogen 

in the galactic plane does not follow the 'Same distribution as the 
atomic hydrogen. Although some caution must be used in using the CO data 

to deduce densities, it seems probable that the H 2 densities are relatively 
larger toward the galactic center than the HI densities and relatively 
smaller toward the outer parts of the Galaxy. Thus, the molecular-to-atomic 
hydrogen ratio would be smaller in the outer Galaxy contrary to the assumption 
of a constant ratio for the local region leading to a decrease in the 7 “^ 3 .y 
emission. Conversely the inverse would be true for the 360° < .^ < 60° region. 

Thus, it seems most reasonable to explain the observed, deviation 
in terms of a combination of a decrease in the molecular-to-atoraic 
hydrogen ratio and in the cosmic ray density in the outer Galaxy and 
a combination of increases in both the quantities in the inner Galaxy, 
rather than variation in either quantity alone. For example, being 
more quantitative, to explain the result on the basis of molecular 
hydrogen alone would require nearly all the molecular hydrogen to vanish . 
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within about 0,5 kpc of the solar system. Since regions B and F are 
generally tangential to the solar circles such a rapid decline in 
molecular hydrogen would imply that the solar system is in a very 
unique position in the Galaxy, a conclusion which is physically un- 
appealing at best. Second, the decline in molecular hydrogen would 
not account for the similar deviation seen in Figure 8 comparing the 
150 MH2 and data. 

The results here then support the concept of the variation of the 
cosmic ray density in association with the matter density on the scale 
of arms and large features, and, as seen in the previous section, argue 
against cosmic ray density variations on a cloud-to-cloud scale, (The 
cosmic ray density is, of course, also believed to vary with height 
above the plane in the manner described previously.) The results here 
do, suggest a larger local matter density than assumed earlier. This 
result would imply a relatively lower cosmic ray density in other parts 
of the Galaxy than used by Kiiffen, Fichtel, and Thompson (1977), since 
the cosmic ray density is normalized to the local matter density. How- 
ever, the increased source function approximately compensates, and the 

general agreement remains * A major contribution from unresolved point 
sources could, of course, affect this discussion. A substantial 7-ray 
intensity from point sources in the inner galaxy would probably weaken 

the argument for a cosmic ray density increase there; however, the con- 
clusion of a low cosmic ray flux in the outer galaxy would seem to remain 
or be strengthened, 

Co The Isotropic Component 

Although there is clearly a substantial apparently isotropic 7-ray 
excess beyond the galactic disk component which was discussed in the 
last section, it should be kept in mind that, with the SAS-2 resolution 
of about 4.5° in the 35-100 JfeV range and the limited statistics. 


isotropy has only been established on a rather coarse scaleo However, 
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the degree of isotropy which has been established, and the steep energy 
spectrum, which also distinguishes it from the galactic disk radiation, 
will be seen to place strong constraints on models for the origin of 
this component. In the next few paragraphs, the galactic and extra- 
galactic hypotheses for the origin of this isotropic component will 
be examined. 

Galactic theories for the origin of the isotropic 7-radiation are 
generally built around a galactic halo which has the shape of a sphere 
or thick oblate spheroid and is populated with cosmic rays held by 
weak magnetic fields. The source of the 7-rays in the halo is Compton 
radiation from cosmic ray electrons interacting with the blackbody 
radiation and to varying degrees, depending on the model, with escaping 
starlight and infrared radiation. The matter and stellar distributions 
appear to preclude bremsstrahlung or point source models, and the 
synchrotron intensity is just too low to contribute significantly to 
the isotropic flxiXo 

Theoretical models Involving a galactic halo have generally 
postulated a halo with dimensions of the order of the galaxy and hence 
a radius, at least in the plane, of about 15 kpc» The implications of 
a much larger halo will be considered later. Since the sun is about 
10 kpc from the galactic center, if such a halo exists and is responsible 
for the 7-rays, a very marked anisotropy would be seen, with the 7-ray 
intensity from the general direction of the galactic axis being much larger 
than that from the same latitudes in the anticenter direction. In fact, 
no such anisotropy is seen; specifically the ratio of the average. 
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intensity in the (300° <^<60°, 20° < jb] <40°) region to that in 
the (100° < & < 250°, 20° < jb] < 40°) region was found, in Section II 
to be 1.10 ^.19 compared to a calculated value. for the same ratio 
of 2,85 for a model with a uniform cosmic ray sphere with a 15 kpc 
radius. - The result for a thick oblate spheroid is similar. . 

It is also of interest to consider the possibility of a halo of 
larger dimensions. If the region is assumed to be spherical and have 
a uniform cosmic ray density, the upper limit (2o) set for the anisotropy 
demands that the radius be at least 45 kpcs. Finally, it is also 
possible to set constraints on the oblateness of such a region which 
might be in the shape of an oblate spheroid. The ratio between the 

absolute value of the polar axis vector and the radius of the circle 
in the plane' is found to be 0.'80jj).15 based on the intensity for 
lb| > 60° compared to that for 20° < jb] < 40° given in Section II. 

This result argues against a major deviation 

from a spherical shape, under the assumptions which have been used here. 

It might be suggested that within the framework of a galactic 
oblate spheroidal or spherical halo, a Compton component might have a 
different and particularly weaker dependence on galactic longitude. 
However, the analysis in Section Illb and related general arguments 
would actually suggest the opposite, if anything, since the cosmic ray 
intensity most likely decreases in the outer regions of the Galaxy and 
increases toward the center. Thus, the Compton disk radiation would 
be expected to show a stronger, not weaker, longitude dependence than 
the conservative (but probably less realistic) assumption of uniformity. 
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This conclusion is probably significantly amplified when diffusion is 
taken into account, .for three reasons; 1 ) an outwardly decreasing 
electron density would result from the diffusion, 2 ) -the diffusion and ■ 
escape parameters would most likely be position and energy dependent 
in such a way as to enhance the anisotropy, and 3) the older electrons, 
which have had the greater energy loss at high .energies tending to 
steepen the energy spectrum and preferentially suppressing the higher 
energy electrons which produce the relevant Compton 7 -rays are more 
likely to be in the outer regions* 

Next, the possibility of a much larger halo can be explored. 

Based on the results given in section II, the electron halo would have 
to have a radius of the order of 45 kpc or more for Compton radiation 
from cosmic ray electrons interacting with the black body radiation to 
explain the isotropic component. Tbts again a conservative limit 
because, as above, it is more likely that the electron density especially 
in the relevant energy range would be larger toward the center; hence, - 
if this possibility is being considered, a substantially larger radius 
would be a better choice for a tentative model* With a dimension of. - 
this size, one is then in effect considering an intergalactic model ■ 
since, if other galaxies are similar, the halos would be beginning to 
overlap. Further since the required density decreases only as the 
inverse of the radius, but the number of electrons to fill the volume 
increases as the cube of the radius, the increase in the required 
number of electrons would be proportional to r . in turn, this demands 
a major increase in source strength or electron age which implies 
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greater energy loss and hence a still larger source® The difficulties 
of the iritergalactic theories (principally intensity) are then rapidly 
approached. 

The differential electron spectrum would have to be quite steep, 
with the spectral index -being- 5+1 for a power law, because of the 
relationship = (m^ +l)/2 between the y-ray and electron spectral 
indices for the Compton radiation. If the escape time from, the halo 
is large compared to 10^ years for most of the electrons in the halo, 
as it probably is since this is a typical estimate for the lifetime 
of cosmic rays in the disk, the spectrum for the electrons with 
energies relevant to 10^ MeV y-rays, that is near 10 MeV, will be 
substantially steeper than that observed in the Galaxy, due to Compton 
energy loss by the electrons. In the case of equilibrium and no escape, 
the spectral index would be only one power greater or about 3=8 (e=g. 
Owens and Jokippi, 1977) however, diffusion considerations, which lead 
to the younger electrons being nearer the disk, combined with the 
likelihood of an energy dependent escape mechanism, make an index in 
the range calculated above, i.e= 4 to 6, seem reasonable. The diffi- 
culties discussed earlier, however, suggest that it is not worth pur- 
suing specific large halo models since only a rather ^ hoc model 
could simultaneously satisfy the intensity and isotropy constraints. 
Among the extragalactic theories, primordial black hole emission 
(Page and Hawking, 1976) universal primordial cosmic ray interactions 
(Stacker, 1971), interactions of cosmic rays leaking from galaxies 
(e=g. Felten, 1973), and emission from normal galaxies (e.g. Strong, 
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Wolfendale, and Worrall, 1976) all face the combined problems of the 
steep energy spectrum and the intensity level. Specifically, if it 
is assumed that normal galaxies are represented by our own, the 
galactic /-ray emission above* 100 MeV would* be inadequate by from over 
one to two orders of magnitude and cosmic rays escaping from the 

galaxies interacting with intergalactic blackbody radiation by an even 
much larger factor. Galactic evolutionary consideration may make normal 

galaxies a somewhat more likely candidate (Licht et alo, 1977); however, 
somewhat more interesting possibilities appear to be radiation from 
exceptional galaxies, such as radio galaxies, Seyferts, and QSO's (how- 
ever, if the radiation originates with electron interactions, the electron 
spectrum must be quite steep), or the y-ray emission (Stecker et al. 1971, 
and Stecker 1977) from the matter-antimatter interactions in the baryon 
symmetric big bang (Harrison 1969, and Omnes 1969). 

The 7-ray spectrum predicted from the latter theory is shown in 
comparison with the data in Figure 11. It must be remembered that, 
whereas the calculation of the y-ray spectrum is straightforward' in 
the baryon- symmetric big 'bang theory, there is substantial question 
as to whether the required separation into regions of matter and 
antimatter which ultimately become superclusters can occur (e.g. 

Steigman, 1974) . In spite of this serious reservation, the very steep 
energy spectrum, which results primarily from the energy shifted tt° 
spectrum from p-p annihilation integrated over cosmological times, 
predicted by this theory Is a strong point in its favor, or at least 
in favor of some cosmological model involving low energy tt'’ decay. 
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The question of radiation from exceptional galaxies is a rather 
speculative one because very little data exist. Until further evidence 
is available, it must remain as a possible explanation; however, the 
combination of the intensity and the energy spectrum remains a matter 
of concern for such theories. Some specific considerations are dis- 
cussed by Worrall (1977) and the reader is, referred to this thesis 
for further amplification of this problem. 

In summary, it seems fair to say that there is no single explana- 
tion for the diffuse radiation which is without difficulties or reserva- 
tions. The status of the experimental data does not seem to justify 
consideration of a combination of several origins in detail for the 
isotropic component, but it is conceivable that both Compton radiation 
from a region of diffusion of cosmic rays from the thick disk discussed 
in sections III -A and III-B and extragalactic radiation of some form 
may be combining to give the observed "isotropic" radiation. 

IV. SUMMARY 

An examination of the intensity, energy spectrum, and spatial 
distribution of the diffuse y-radiation observed by SAS-2 in the energy 
range above 35 MeV has shown- that it consists of two components. One 
component is generally correlated with galactic latitude, the atomic 
hydrogen column density, and the continuum radio emission. It has an 
energy spectrum similar to that in the plane, and joins smoothly to the, 
intense radiation from the plane, which has been examined in some detail' 
previously. It is therefore presumed to be of galactic origin. The other 
component is apparently isotropic, at least on a coarse scale, and 
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has a quite steep energy spectrum. The features -of the galactic 
component of the y-radiation and conclusions which“can he inferred 
from the data are now summarized first: 

•There is a strong linear correlation between both the 35 MeV 
< E < loo MeV and the > 100 MeV y-ray intensity and the radio 
continuum radiation as measured at 150 MH^. 

•For |b|^12,8°, there is a good linear correlation between the 
y-ray intensity in the same two energy ranges and the deduced 
atomic hydrogen column density. The constant residual y-ray' 
intensity at the extrapolation to zero column density implies 
an "isotropic” y-ray component, whose properties are summarized 
below. 

•For jb| < 12.8° and 60° < I < 300°, there remains a correlation 

between the y-ray intensity and the atomic hydrogen column density 
deduced from the 21 cm emission, but there is less y-radiation 
than would be expected from the linear extrapolation based on 
data for jbj > 12.8°. The magnitude of the decrease implies -that 
most likely both -the cosmic ray density and the molecular -to -atomic 
hydrogen ratio decrease relative to the local region for galactic 
radii larger than that for the solar system, although the 
possibility that only one of these two components decreases 
dramatically cannot be excluded. 

•For 300° < I < 60°, the y-ray intensity from the galactic plane 
actually is well above that suggested by a linear extrapolation 
based on data for|bj > 12.8°. This result implies that either 
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the cosmic ray density and. the molecular— to-atomic hydrogen, 
ratio both increase in the principal features in the inner 
galaxy or one of the two Increases dramatically. Previous 
studies of the galactic plane along suggest that the most 
plausible assumption is that both quantities increase. 

•The energy spectra of the ytnys from the galactic plane and 
from the high latitudes are the same within uncertainties. This 
similarity and the smooth connection in intensity suggest a 
similar origin. 

•There is a larger intensity in the 35 MeV to 100 MeV range 
relative to that above 100 MeV than expected from most current 
estimates of the cosmic ray electron spectrum. If this 35 MeV 
to 100 MeV excess emission is interpreted as being due largely 
to cosmic ray electron bremsstrahlung, a significant C'^J3 times) 
enhancement in the low energy (35 MeV to several hundred MeV) 
cosmic ray electron spectrum over previously assumed spectra is 
implied. If the excess is attributed largely to Compton 
radiation, the electron intensity above 50 GeV must be sub- 
stantially above that given by current measurements at the earth. 
No enhancement is necessary for the 1 to 30 GeV range, so 
there is not necessarily a conflict with the radio continuum 
observations. A much larger low rigidity solar modulation would, 
however, be implied for the cosmic ray electrons. 

The "isotropic" component was found to have the following 
properties and implications: 
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•It appears to exist for all seven regions of the sky examined". 

• -5 -2 

•The intensity above 35 MeV is (6.3 + 1.4) x 10 photons cm ' 
ster ^ s ■ 


•The energy sp.ectrum is quite steep, having a 2,.:85 + 


0.50 

6.35 


differential power law index for a best fit straight line 


between 35 MeV and vl50 MeV. 


•When extrapolafed to 10 MeV, this deduced power law joins 
smoothly to the diffuse "isotropic" intensity- measured" at 
low energies (0.1 to 10 MeV) . 

•Relatively few theoretical models are consistent with the 
combination of the observed spectral shape intensity and lack 
of strong anisotropy, and those are not entirely free of other 
difficulties o 


°The origin of this "isotropic" y-ray emission’ must, however, 
be considered a major open question with further high sensi- 
tivity measurements of the Intensity, uniformity, and energy 
spectr"um being very important o 



39 


Table I 

> 

Comparison of Predicted* and Observed Gamma Ray Intensities 


’20 2 

A. Polar Region (typical) = 2,4x10 atoms/cm ) 


Energy Range 

(35-100 MeV) 

(>100 MeV) 

Cosmic Ray Nucleon, Matter Interactions 

(0.08-0. 12) xl0“^ 

( 0.32-0. 48) xlO"'^ 

Cosmic Ray Electron, Matter Interaction 

(0.l4-0.33)xl0“^ 

( 0.08-0. 20) xl0“^ 

All Compton Interactions 

(0.08-0, 32) xl0“'^ 

(0.05-0.20)xlb“^ 

Total Predicted 

(0.30-0.77)xl0“^ 

( 0.45-0, 88) xlO"^ 

Observed"*" 

(0.83+0.21)xl0“^ 

(0.76+0.l9)xl0“'^ 

• » « 

• •» 

B. |b| = 15° typical) = 11x10 

^*^atoms/ cm^) 


Energy Range 

• (35-100 MeV) 

(>100 MeV) 

Cosmic Ray Nucleon, Matter Interaction 

(b.37-0.55)xL0“^ 

(1.47-2. 20) xl0~^ 

Cosmic Ray Electron, Matter Interaction 

(0,64-1. 51) xl0“^ 

(0.20-0.80x10“^ 

All Compton Interactions 

(0.31-1.24)xl0“^ 

(2.0-3.9)xl0“'^ 

Total Predicted 

(1.3-3.3)xl0“^ 

(2.0-3.9)xl0“^ 

Observed^ 

(4.0+0.8)xl0"'^ 

(3.6+O.'7)xl0“^ 


*A range' of predicted values is given, rather than a specific 
value, as discussed in the text. 

“ 1 “ 

The errors quoted for the observed intensities reflect systematic 
uncertainties and the estimated uncertainty in separating the two 
components, which for example, do not effect ratios of total Y“tay 

.intensities in the same measured energy intervals . 
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FIGURE CAPTIONS '■ ■ 

Figure 1 Diagram of the sky in galactic coordinates showing the. 

portion of the sky viewed by SAS-2 (shaded) and the 
regions -used- in- -the present -analysis-, which- are -bounded- by 
dark lines and designated by letters. 

Figure 2 Distribution of v-ray (E ^ 100 MeV) intensities as a 

function of galactic latitude for the indicated longitude 
intervals. Points for |b! > 61.0° represent data summed 
over all longitudes and are Indicated by open circles. 
Latitude division boundaries are not necessarily Integers 
because the regions represent sums of smaller areas formed 
from dividing the sky into equal solid-angle regions formed 
by fixed longitude intervals and 72 latitude intervals. The 
intensities shown are based on assuming a typical spectrum; 
however, since the actual values are somewhat affected by the 
energy spectral shape because of the finite energy resolution, 
the reader is referred to the later parts of this paper for 
improved estimates of the intensity and energy spectrum based 
on a two-component analysis and specific energy spectra. 

Figure 3 Distribution of y~ra.y (35 MeV < E < 100 MeV) intensities 
as a function of galactic latitude for the indicated 
longitude intervals. Points for |b| > 61,0° represent 
data summed over all longitudes and are indicated b}' open 
circles. See also the comments in the caption for Fig. 2 
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Figure 4 


Figure 5 


Figure 6 


Distribution of y-ray (E^ ^ 100 MeV) intensity as a function 

of atomic hydrogen column density deduced from 21 cm radio 

data for |h | > 12.8*^ for the indicated longitude and 

latitude intervals. Error bars shown on the y-ray Intensities 

are statistical. An uncertainty of approximately 10% should 

be attached to the 21 cm column densities. The_ solid, line 

is the best fit to all the data points shown. Upper right: 

2 

acceptable (Icr above the minimum 'f, ) values of the fitted 

parameters in the equation + B. 

Distribution of y-ray (35 MeV < < 100 MeV) intensity as 

a function of atomic hydrogen column density deduced from 
21 cm radio data for |b| >12.8° for the indicated longitude 

s 

and latitude intervals. Error bars shown on the y-ray 

intensities are statistical. An uncertainty of approximately 

10% should be attached to the 21 cm column densities. The 

solid line is the best fit to all the data points except the 

point at (1.8, 1.1) in the region 250° <Z< 300° (this 

anomalous point is incompatible with any linear fit to the 

rest of the data). Upper right: acceptable (la above the 
2 

minimum f. ) values of the fitted parameters in the equation 


I = A N,,,. + B. 
y HI 

Distribution of y-ray 100 MeV) intensity as a function 

of atomic hydrogen column density deduced from 21 cm radio 
data for jbj <60 for the indicated longitude intervals. The 
longitude intervals not shown did not have meaningful data 
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for the present purpose because of knovm ooint’ source 
contributions. The’ solid line is the one derived from the 
high'-latitude data shown in Fig. '4. 

Figure 7 Distribution of y-ray (35 MeV< < 100' MeV) intensity as 
a function 'of atomic hydrogen column density deduced from 
21 cm radio data for [b| < 60° for the ind'icated longitude 
intervals. The longitude intervals not shown did not have 
meaningful data for the present purpose because of known 
point source contributions. The solid line is the one 
derived from the high-latitude data shown in figure 5. 

Figure 8 150 MHz brightness temperature as a function of atomic 

hydrogen column density for |bj <60° and the indicated 
longitude ranges. The estimated uncertainties are approximately 
10% for the 21 cm column densities and 15% for the 130 MHz 
brightness temperatures. The solid line is the best fit 
to all the data with jb| >12.8°, averaged over the same 
regions used for the y-ray studies. Inset; acceptable 
(la above the minimum ) values for the fitted parameters 
in the equation = A + B, using only data with 

|b|> 12.8°. 


Figure 9 Gamma-ray intensity (35 MeV< E< 100 MeV) as a function of 

150 MHz brightness temperature. Only typical error bars are 
shown because of the high data density of the figure. The 
estimated uncertainty inthe 150 MHz brightness temperatures 
is approximately 15%. Straight line is the best fit to all 
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the data points. Inset: acceptable (lo above the minimum 

2 

) values of the fitted parameters in the equation I = 


^ ^150 


Figure 10 (35 MeV< E< 100 MeV) y-ray intensity as a function of (E> 100 MeV) 
Y~ray intensity. Only typical uncertainty limits are shown 
because of the high data density of the figure.- The high 
energy y-ray data have been ’’defocused" for comparison to the 
35-100 MeV data, as described in the text. The straight line is 
the best fit to all the data points. Inset; acceptable 
(la above the minimum Si ) values for the fitted parameters 
in the equation ^ +B. 


Figure 11 Gamma-ray spectrum for the ’’isotropic” component of this work, 
together with other data and the theoretical spectral shape 
predicted by the baryon symmetric big bang theory discussed 


in the text. 
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